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; Abstract 

(N : 
> ■ 

' It is known that gauge fields defined on manifolds with spatial boundaries 

^ ■ support states localized at the boundaries. In this paper, we demonstrate 

. how coarse-graining over these states can lead to an entanglement entropy. 

\ In particular, we show that the entanglement entropy of the ground state for 

' the quantum Hall effect on a disk exhibits an approximate "area " law. 

Oh: 

(D . I. INTRODUCTION 

; 

\ Recently there has been renewed interest in the origin of the so-called area law for 

^ ■ black hole entropy. One reason for this interest is that it can shed some light on the possible 

■ substructure for a quantum theory of gravity. In this regard one often invokes the hj^othesis 
that there are excitations of the black hole horizon leading to surface states associated to 
the horizon. In a previous paper using a very simple treatment of the constraints of 
canonical gravity, we have shown that these states arise in a manner similar to the edge 
states in quantum Hall effect (QHE) IQ-Q. We also speculated there that tracing over these 
edge states might lead to an entropy for black holes as happens with the calculation of 
entanglement entropy |^,|^. 

A complete calculation to verify the above speculation is a rather formidable task. How- 
ever, as a step forward, we can try to find evidence in other models where edge states arise 
and lead to important physics. With this in mind, in this paper we study the effective field 
theory describing QHE 0] on a disk and calculate the entanglement entropy [^,^ that arises 
when we trace out the edge degrees of freedom. More specifically, we see that the ground 
state of the system turns out to have non-trivial correlations between the bulk and the edge 
degrees of freedom. For this reason, tracing over the edge degrees of freedom leads to an im- 
pure density matrix. The entanglement entropy that we calculate is the entropy associated 
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with this density matrix. The result we find is quite interesting - the entaglement entropy 
obeys an approximate area law [ that is the entropy is approximately proportional to the 
perimeter of the disk ] for weak coupling between the edge and the bulk. An area law of 
the similar type has also been discovered for the thermodynamic entropy (2+1) dimensional 
black holes recently 0,0]. 

The paper is arranged as follows. In Section 2, to keep the paper self-contained, we discuss 
how edge states arise for gauge theories defined on manifolds with spatial boundaries. In 
Section 3, we discuss our model , namely Maxwell-Chern-Simons ( MCS) theory defined on a 
disk coupled to a chiral scalar field living on the edge of the disk. This is a phenomenological 
model that describes QHE (along with the edge currents present in a Hall system). In Section 
4, we find the Hamiltonian for the system by eliminating the unphysical modes using the 
Gauss law. This Hamiltonian is quadratic and hence in principle, one can find ground state 
wavefunction [} In Section 5, we calculate the ground state using perturbation theory and 
the entanglement entropy associated with this ground state. We find that, in the above 
weak-coupling limit, this entropy is approximately proportional to the perimeter of the disk 



( "area law"). This result can also be checked using the so-called replica trick |]T0|, which 
we discuss in the appendix. 



II. EDGE OBSERVABLES 

In this Section, we will see how edge observables occur for gauge theories on manifolds 
with boundaries. We will take the MCS theory on a disk as a typical example since we 
will be using this model also in the later Sections. 

Here, we use the following conventions: 

1. Greek and Latin indices take values 0,1,2 and 1,2 respectively. 

2. The three-dimensional metric r/^i, is specified by its nonvanishing entries 

r]oQ = —1, 7]ii = ?722 = +1 while the three-dimensional Levi-Civita symbol is e^'^^ with 
e0i2 = +1. 

3. The spatial metric is given by the Kronecker 6ij symbol while the two-dimensional 
spacetime and spatial Levi-Civita symbols are e^*^ and e*-' = e"*-'. 

We will also assume that the disk has a circular boundary and radius R. 
The MCS Lagrangian is 

L = f SxC 
Jd 

C = -\ F,,F^- - ^ t^'^A^d^Ax . (2.1) 
The Poisson brackets (PB's) for ( p.l|) are 



We thank R. Sorkin for emphasizing this point to us. 
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{A(x),A,(t/)} = {n,(x),n,(i/)} = o, 

{A,{x),Uj{y)} = 5,,5\x-y). (2.2) 

Here and in what follows, all fields are to be evaluated at some fixed time while the PB's are 
at equal times. Thus x° = in (|2.2| ). As usual these PB's will be replaced by commutation 
relations (CR's) in the quantized theory. Also A and 11 give the magnetic field B and the 
components := Ei of the electric field by the formulae 



B — eijdiAj, 

-(n. + ^ 
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Ei = -(n, + -^ei.A,). (2.3) 



The Hamiltonian and Gauss law for ( |2.1| ) in quantum theory are 



H = [ dPxU, 

JD 
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G(x)|-)=Oforx|aD = 0, (2.4) 

where 

G{x) = -J^ d^x da^'^ [n, - ^e,,A,] (2.5) 

and |-) is any physical state. 

The reason for the boundary condition on x is that it is only with this condition that 
vanishing of G{x)\') is implied by the usual Gauss law diijli — ^tijAj) ^ (after a partial 
integration). Note that it is necessary to rewrite this classical Gauss law in terms of G{x) in 
quantum theory. This is because Ilj — ^eijAj are operator valued distributions in quantum 
theory so that their derivatives have to be interpreted by smearing them with the derivatives 



of suitable test functions |]T2[. But it is only when the test functions satisfy the above 
boundary conditions that the smeared Gauss law restricted to the classical context vanishes 
by the classical Gauss law [0 . 

The edge observables Q{A) are obtained from G by changing the boundary conditions 
on X- They are 

Q{A) = — J (fxdiA[Ili — ^^eijAj] , AI^d = not necessarily . (2.6) 

They are the generators of the affine Lie groups LU{1) and have the commutators 

[Q(A),Q(A')] = -lan [ d^xeijdiAdjA' = -lan I AdA'. (2.7) 

The action of Q{A) on |-) depends only on the boundary value of A since the difference 
between Q{A) and Q{A') when A and A' coincide at the boundary is a constraint and hence 
annihilates the physical states. Furthermore, it commutes with observables localized within 
D, its action on these observables being that of a G{x) Hill- Hence it can be regarded as 
localized at the edge. 
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III. MCS ACTION WITH THE CHIRAL SCALAR FIELD 



The MCS action we have considered till now suffers however from an "anomaly" . Under 
the gauge transformations 

Ai + diA, 

^ E,, (3.1) 
for arbitrary A, the MCS action Sbuik is not gauge invariant: 



Sbuik ^ Sbuik - / d'xe^'d^AA, (3.2) 

2 JdM 

{Ai being D x H). If we require that physics be gauge invariant, then we must modify the 
MCS action Sbuik by a surface term as follows: 

Stot = Sbuik + ^ / d'^xe'^'d^M. I rf2x(D^0)(D», (3.3) 

2q JdM oTT JdM 

D^<P = d^<j) - qA^. (3.4) 

Here is a positive constant and q is the charge by which the field living at the boundary 
is gauged. Also d'^x = dt Rd6, R being the radius of dD and 6 G [0, 27r] its angular variable. 
Under the gauge transformations ( p.l|) , transforms as 

0^0 + gA, (3.5) 

so that $ = e**^ transforms like a complex scalar field with charge q. [Since $ is the true 
charged excitation, we also have the identification 0^0 + 27r.] 
With this transformation law for 0, we see that 

Stot ~^ Stot- (3-6) 

So far, the coefficient in ( p.3|) is arbitrary. However, if this action is to describe QHE 
T3|,|3[|, then the currents due to the edge scalar field are required to be "chiral". In this 



case, it can be shown ^ that the coefficient gets uniquely fixed according to 

f = 2vr^ (3.7) 
q^ 

With P fixed in this way, we can impose the "chirality" constraint that the edge fields 
are left-moving, namely that 

Do(f) - Decf) = 0. (3.8) 
Thus the total action that correctly describes QHE is 

Stot = Sbuik + ^ / d'xe^-'d^Mu - ^ / d'x{D,4>){D^'4>). (3.9) 
Zq JdM 'iq JdM 
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IV. THE HAMILTONIAN AND QUANTIZATION 



In this Section, we will find the Hamiltonian for the action (|3.9| ) and then set up the 
formalism for the quantization of this interacting system. 

From the action (|3.9|) , we find the momentum conjugate to Ai to be 

n, = tEi - -feijAj (4.1) 
(where Ei := Foi) and the momentum conjugate to (f) to be 



The Hamiltonian for this system is 

Htot = Hbulk + Hedge 

1 /■ r2g%T-r , ^2 , (^H 



and the Gauss law constraint is 

<3{X) ■■=- I c/\9,x(n, - ^e,, A,) - g / i?c/^ x(n^ + ^00 ^ 0, (4.4) 

JD JdD /Q' 

R being the radius of the disk and ^(mod 2tx) the angle coordinate for dD. To quantize the 
system described by ( |4.3|) and ([4.4|) , we now need to choose mode expansions for the fields 
Ai, Ilj, and H^. As usual, the coefficients of these modes will play the role of creation and 
annihilation operators in the quantized theory. 

In an earlier work ||TT],|T2[, we have quantized the bulk theory alone (namely the MCS 
action without the chiral scalar field). 

There we showed that among the various permissible mode expansions parametrised 
by a real parameter A, only the one characterized by A = allows the existence of edge 
observables. This is because each A specifies a domain Vx for the space of one-forms [Ai or 
Ilj) and the edge observables do not leave this domain invariant unless A = 0. Since edge 
observables are crucial in our work here, we will use only the mode expansions specified by 
the choice A = 0. 



The corresponding mode expansions for A, 11 and E are 



A = [anm^T + '^nrn 7^ + an^=)+ h.C, (4.5) 



n = {jlnm^t^nm + TT^lVt * + Vn^OuK) + h.C, (4.6) 



^ = '-^""■^ ^ ^"'""V?'' ^ ^ ^«n)^n + h.C. , (4.7) 
_ ^ {*) 

p\ > = TT^ ' A- n 
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summation over repeated indices [ n over positive and m over non-negative integers] being 
understood. Here, we have used ( |4.1|) and also the notation of forms while * refers to 
the Hodge dual [M. Further we have scaled the modes of lO] for later convenience and 



dropped the superscript (1) employed in |jTTI] to emphasize one-forms. The modes used above 
are defined as follows ( see for details): 

K ■■= cifrV"^). (4.8) 

The d refers to the exterior derivative, the J„'s refer as usual to the cylindrical Bessel 
functions [|T5| while the Wnm^s are such that Jn{uJnmR) = 0. Also the A^nm's are normalization 
constants chosen such that ^nm * ^nm = — 1, bar denoting complex conjugation. 

In ( |4.5| -C7), it is understood that n is summed over all non- negative integers for the first 



two terms and all the positive integers for the third term while m is summed over all positive 
integers. 

From the PB relations ( p.2|) and using also |2.3|) , we get the following non- vanishing CR's: 



r (*) (*)t 1 _ r (*)t „(*) 1 _ r„(*) „{*)t i _ r„(*)t J*) l _ ,^ , 

fp pWt 1 _ U J*) 1 _ ,X , 

= ["LPn'] = iSnn'- (4.9) 

We also mode expand the edge fields cj) and H,^ as follows: 

^-^^ ^ e*"^)+ h.c, (4.10) 



H, = (vr„^^?^e^"^) + h.c. (4.11) 



We have suppressed the winding modes in writing ( [4.10 , 4.11 ) as they are not important for 



this paper. Also n here is summed only over non-negative integers. Though ( [4.10D looks 
singular for n = 0, it is all right for our purposes since what occurs in the Hamiltonian and 
the Gauss law is 0' (and not 0) and this itself admits a well-defined mode expansion even 
for n = 0. Once again, the only non-zero CR's are 

[0„,4,] =z5„„/. (4.12) 

The Gauss law ( [4.4] ) in terms of these modes takes the form 

(e« - ^an„.)|-) = (eW - ^a_)t|.) = (4.13) 



and 



[{Pn - '-a^) + (vr„ + '-<Pn)] I-) = [{pi + '-ai) + (vrt - l^t )| |.) = 0, (4.14) 
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|-) being any physical state. Here ([4.13|) is the condition on the modes *'^nm while ( [4 .141) 
arises as the constraint on the modes hn,hn- 



We can use ( [4.13| ) to eliminate e*^ in favour of a„m- In that case the Hamiltonian ( |4.3| ) 



in terms of these modes becomes 

^edge 



Htot = 


Hhulk + He 


Hhulk = 


l|et e 


Hedge 


n i 

(7r„ + -(pn 


where 




and 





nm 



f^)2N t Q \ , ^fc^C +C C^l 

V nt 



^^^'^-t^nm'A^nm'-C(^nm'-R)«nm') + (^Tn - ^0n)^(7rn - ^0n)] (4.15) 



Cn = (Pn + -an)^ (4.16) 



4(^) = ^Ux). (4.17) 

If we are in one chiral sector of the edge theory (the physically relevant sector when dealing 
with QHE), then we can also impose the condition 

(^n-^0n)|-) = O (4.18) 

on physical states |-) (such a condition arising from the mode expansion of the chirality 
constraint (^.8[)). 

We can now also use the Gauss law ( [4.14| ) to eliminate (7r„ + |0„) in terms of (p„ — |a„). 



The edge Hamiltonian then reduces (in the above chiral sector) to 



Hedge — ^[(^'" + 2*^" ~ ^\ ^nniNnmJ'n {^nmR) O-nm) ^ 



[Pn + ^"n - -Ry ^^^^nm'A^nm'-C(^nm'-R)anm')]- (4.19) 

Therefore the task of quantization reduces to the diagonalization of the system described 
by the following Hamiltonian : 

H=-[e^ e + {—?)a) a \ + —c^c + 



n , i 

;^(4-^ 



]j ^nmNnmJ'n{^nmR)a'nm){Cn — R]j ^^ -Unm' Nnm' J'n{^nm' R)(^lirn') ■ (4.20) 



The Hamiltonian in ( [4.20|) has also been normal ordered, the vacuum |0) being defined by 
enm|0) = OnmlO) = c„|0) = 0. Notc that H preserves the constraints ( |4.13| ),( ^?T^ ) and ( [4.18[ ). 



V. THE GROUND STATE AND ITS ENTANGLEMENT ENTROPY 



We are not interested here in finding the entire spectrum of the Hamiltonian (|4.20| ). Our 
objective is in finding the ground state of the system. As is clear even from ( [4.20| ), this 
system has correlations between the edge degrees of freedom a„ and bulk degrees of freedom 
Cnm, cinm- The density matrix obtained by tracing out the edge states will therefore be 
impure and will have a non-zero entropy associated to it. 

We will now find the ground state of this interacting system by using perturbation theory 
with q being the perturbing parameter. Before starting the calculations, it is important to 
note that in physical situations in QHE, the coefficient an can be written as kq"^ where k is 
a number of order unity. We will therefore replace an in ( |4.20|) by kq"^ to make the order in 
q explicit: 



An \ ^ r 4- 



W t 



n 
R 



fO/Trk fO/Kk 

~ ^nm\/^nmi the 



Keeping terms only to order and defining e„m = e„m/ yj^nm and a 
above Hamiltonian can be rewritten as 



,27ikR 



t 



We now define annihilation operators Anm and Bnm as follows: 



(5 



Br. 



1 

V2 



Their non- vanishing commutators are contained in 
In terms of these variables the Hamiltonian is 

X Ti 

H = —OJnml'Anm'^nm, + Bl^^Snrri] + ^'"n''"- 



'/ Tinkuj 



"-^nm<-^n(^nm-R) (Cn(^n77j •^nm) ~l~ '-'n('^nm ^nm)) ~l~ 



(5.3) 



(5.4) 



i V t 

2 k irkR j 

?^[jc]jC„ H — y/ UJnm'^nm' NnmNnm' J'n i^nmR) J'n i^nm' R) {-Anm-^nm' + ^nm'^" 



(5.5) 
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upto additive constants caused by ordering of operators. 

We can now find the ground state \'^)gnd to order g^. For convenience let us fix the 
angular momentum n. The true ground state will be a product of wavefunctions with 
different n's. Therefore, 

= ni*'^)fl"'i' (5-6) 
n 

where, 

\ |n\ri ^ 7^nkuJnn,N^m{Jn{^nmR)y ^ , 



|1)cJ1)^„^WA/ 7 n I ) + 

m V t 2 



2 / R '^nm 



2 |-|\ |.|\ .Tmk,\/ 2(jJn^Unm' ^nm^nm' Jni'^nmR) Jni ^nm' R) 

m,m' Vi?"r2/Vii"r 2 / 

Here |0) refers to the unperturbed ground state : 

Anm\0) = Bnm\0) = CnlO) = 0. (5.8) 

Also the state refers to the first excited state created by the creation operator Al^^, 

a similar notation being used for the other states. 

Prom this ground state, we can also now define the corresponding density matrix p and 
then obtain the reduced density matrix by tracing over the bulk states. [Actually, we may 
want to trace over the edge states if the objective is to calculate the entropy due to a lack 
of knowledge of the edge states. But the answer for the entropy itself is the same whether 
we trace over one or the other.] To order q^, we have for this reduced density matrix, 

Pred := Tr^bulk){p) = |0)a„a„(0|(l - P„) + |l)„„a„(l|p„ (5.9) 

where 

m ^ yR^ 2 ) 

The dependence on R factors out in the expression for p„ as can be seen by rewriting it as 
follows: 

p„ = , i-nk)jY. (n+^)^ 

~ '''' ^ tr (2n + Xnm)2 

^ q 7 2. , ^ ^2^ large Xnm := ^nmR, (5.11) 



9 



Xnm being zeroes of J„. This shows that p„ is proportional to R. Note that pn is finite. 
This is because for large m, 



\Xr, 



TT 1 

~ -(2m + n- -), 



NnmJ'niXnm) ~ 0{ ) 

"X.nm 



SO that 



(2n + Xnmf 



(5.12) 



(5.13) 



Here the above behavior of Xnm, and the behavior of JniXnm) and JniXnm) for large \xr. 
are known from the literature ||16| to be 



-^cos(x„™-^-^) = 0( ^ 

" A.nm ^ ^ 



2 mx TX 

sin(x„m - ^ - -^)=^ 



\/Xnm 
, 1 



2 4' 



while that of AL™ is obtained from the normalization condition 



D 



Xnm ji2 { XnmT \ , ^ j2 { XnmT 



R "^""^ R 



R 



rdrdO = 1 



on \l/nm- Hence, 



27r 



du 



,n 



+ {-Jn{u)y 

u 



u du. 



where u = It shows that depends only on Xnm and gives 



9{w- 

-* ^ nr, 

dXnr 



constant for large \xr. 



Hence 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



Ar_, = 0( 



as \Xn 



oo 



Putting these estimates together we get (|5.13|) 

We can now infer that p„ being finite, the entropy of the rath mode 

Sn ■= -Pn logPn - (l - P„.) log(l - p„) 



(5.18) 



(5.19) 



is also finite for generic p„. 

[However, the total entropy defined as S* := J2n Sn can diverge unless one imposes a 
cut-off (for example, on the maximum allowed n).] 
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VI. CONCLUDING REMARKS 



In this concluding section, we draw attention to the approximate "area law" obeyed by 
the above entanglement entropy. From the expression (|5.19|) for the entropy Sn, we see that 



in the limit of q being small, we can approximate it by just the first term. Namely, 

Sn^Pn- Pn log p„ for q Small (6.1) 

Since we know from (|5.11| ) that p„ itself scales as R, we see that S'„ also scales as R (apart 
from the logarithmic term). For a disk, this is therefore a statement that the entropy 
approximately scales like the perimeter of the disk (or the area of the boundary of the disk) . 

It is interesting that such an "area law" has emerged for the entanglement entropy of 
edge states. Before we draw general conclusions, however, we should check whether such a 
law is a mere coincidence for the 2+1 theory describing Hall effect or whether it holds more 
generally. One possibility is to find the entanglement entropy for a similar theory in 3+1 
dimensions, but without a Chern-Simons term. If the "area law" does hold in such more 
general situations, it seems reasonable to conclude that the black hole entropy is at least 
partially the creation of the entanglement entropy arising out of a lack of knowledge of the 
edge states describing excitations of the black hole horizon. 
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APPENDIX: THE REPLICA TRICK IN ACTION 



In this section we show that the results given above also follows from another method, 
namely the replica trick |Ty]. Our Hamiltonian, given by (|4.2CI| ), is quadratic in the various 
field modes and their conjugate momenta. The ground state wavefunctional for this system 
will then be a generalized Gaussian. To simplify the problem and to highlight the important 
aspects, we will restrict ourselves to one tower of radial modes with any fixed m, which 
we will call a„, and keep all the edge modes a„. The measure defining the scalar product is 
Yin dandaldandaj^. The calculation is performed using the Schrodinger representation where 
we have the substitutions 

■ 9 . . d , , ^, 

Pn -i—j , Pk -^7^- (A2) 
dal oan 

The ground state wavefunctional is found by using the following Gaussian ansatz which is 
consistent with angular momentum conservation: 

*({a„}, {4}; {an}, = J] Ar„e-(^"""'^^+^"""""+^"'^"'^"+^"""'^"). (A3) 

n 

Here Hn-Mi is the normalization factor for the wavefunction and An, Bn, Cn, Dn are real 
constants to be determined. It is also understood that f denotes complex conjugation here. 
The Schrodinger equation is 

H^{{an}, {an}) = E^{{an}, {an}) (A4) 

where our Hamiltonian is given by ( [4.20[ ). The constants An, Bn, Cn and Dn in the wavefunc- 
tional can be determined from { \K^ ) and requiring the wavefunction to be normalizable. 
An, Cn and Dn turn out to be 



C 



l2naH 2inUJnNnmJniX' 



nm I 



V nt {2En + ^ + ^) ' 
Dn = 0, (A5) 

while Bn is determined by a complicated cubic equation. The solutions to this equation are 
rather cumbersome and are not illuminating. However, in the weak coupling — 0) and 
large radius limit ( R oo), one gets the unperturbed value, 

Bn = OOn- (A6) 

The ground state energy is given hy E = En where 

En={^ + ^)An + \Bn, (A7) 
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Hereafter we will write our wavefunctional as 

*({a„}, {4}; {an}, {4}) = J] Ar„e-(^"""""+^"""""+^"'^'^""\ (A8) 

n 

as Dn = 0. 

Given the above form for the wavefunctional, one can readily form the ground state 
density matrix p : 

p({a„}; {an}\{aj; {a'J) = *({a„}; {«„})**({<}; {a'J) 

— ■J~j' |_yy'^|2g-{A„(a„a3i+a^a^t)+-Bn(a„a3i+aJja^^)+C'„(4iQ!„+Q^^aJi)} 
n 

= l[p^''\an,anK,a'^). (A9) 

n 

The reduced density matrix is found by tracing over the set of variables which is not 
being observed , say a„'s. Therefore, the reduced density matrix in this case would be given 
by 

Pred.{{an}] Wn}) = IT / danddnp^"-\an] an\an] a'J. (AlO) 
A straightforward Gaussian integration over the complex variables a„ leads to 

n 

= pS(an; a'J (All) 



where Af'n = -^n^J-^ and 7^ = -a^- The normalization factor N' can be computed easily 
by recalling that for any density matrix p, we require tr p = 1. This gives 

\M'n? = -{2An-^n) (A12) 

The entanglement entropy S = — tr pred In Pred can be computed from the reduced density 
matrix using the so-called replica trick according to which 

S = -^{tr {pZd))\N=u (A13) 



Pred here being the operator with the kernel ( [A10| ). 



Now, one can find for real and positive A the following equation 



Prediian}; {a'n}) =11/11 dZi^ndZi^npi^li^^n, Z^^n) P^^edi^l^n, Z2,n) " " " pS(^W-l,n; "n) 
n •' i=l 
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n 

. . . X y dZl dZl e~2"^"^l."^l."+T"(^l."""+^2,n2:i,n) 

(which on integrating over the 2:1 „ variables becomes ) 

where we have repeated used the identity 

dzdze-^'^+^'+''' = ^e^. (A15) 

The bars here denote complex conjugation ( just as the daggers). 
Hence, we can see 



(Pred) = n / dandalp^J^^aian, an) 

n 



Tc . 1 vr 



^94 ^ 94 — 'v f Tn \N-i 
(OA --v 1^ fl - i^yl)^ 



where we have used ( |A12| ). This expression is invariant under the simultaneous rescaling 
Bn, Cn — > AA„, A-B„, AC„ of the constants. Hence so will be the entropy. Using it in 
( [A13|) , it readily follows that 



^ (2A„ln2A„-7,ln7„) -ln(2A„-7„) 



(2A„ - 7„) 



(A17) 



The constants Cn are determined by the coupling between the variables {a„} and {ctn}- 
When the couplings vanish, C„ are zero and so are 7„ and hence, as one can see from ([A17D , 
the entanglement entropy vanishes as well. 

In the weak coupling limit (g' — 0), the following simplifications occur. One can write 
using equations ( |A7|JA5| ) 



Cn 



nn 

~t {ujn + '-^ + f ; 



/27rA;g2 2inuJnNnmJnix 
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where we have dropped the term proportional to in the denominator as it would be small 
compared to the other terms in the weak coupling limit. 
Therefore, one finds 



In 



nmjXn 



2B„ 



iXn + 2ny 



(A19) 



where we have used the relation Xn = oJnR- Note that 7„ is nothing but p„ for a single bulk 
mode. Due to the relation An = ^, the entropy expression now reduces to 



(1 - In) 



-7„ln7„) - ln(l - 7„,) 



(A20) 



For small 7„ this expression reduces to that appearing in equation ( |6.1|) and hence shows 
exactly the same scaling behavior. 
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